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We consider several observers who monitor different parts of tlie environment of a single quantum 

system and use their data to deduce its state. We derive a set of conditional stochastic master 

equations that describe the evolution of the density matrices each observer ascribes to the system 

under the Markov approximation, and show that this problem can be reduced to the case of a single 

"super-observer" , who has access to all the acquired data. The key problem - consistency of the sets 

of data acquired by different observers - is then reduced to the probability that a given combination of 

data sets will be ever detected by the "super-observer" . The resulting conditional master equations 

^^ , are applied to several physical examples: homodyne detection of phonons in quantum Brownian 

f^ • motion, photo-detection and homodyne detection of resonance fluorescence from a two-level atom. 

(^ ' We introduce relative purity to quantify the correlations between the information about the system 

^S| , gathered by different observers from their measurements of the environment. We find that observers 

gain the most information about the state of the system and they agree the most about it when they 

measure the environment observables with eigenstates most closely correlated with the optimally 

predictable pointer basis of the system. 
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^^ Information about the state of a quantum system is usually obtained not from direct measurements on the system, 



^ ■ I. INTRODUCTION 

> 

^\} ■ Information about tfie state of a quantum systei 

^— V , but rather by monitoring its environment [l|,y,|3{. Therefore, the environment is not only a reservoir that selectively 

r — ■ destroys quantum coherence, but also a "communication channel" through which observers find out about the system 

f^ . Q- The formalism that ascribes to the system a time-dependent state deduced from a complete measurement of 

the environment - a quantum trajectory - has been introduced some time ago IM la S HI- Persistent monitoring 

of a quantum system by the environment can single out a preferred set of states, known as pointer states, which 

are the most robust in spite of the interaction with the environment, that is, least perturbed by it 0,[i3|- In uM 

^ i' we showed that under reasonable assumptions pointer states remain the most robust even when a single observer is 

performing continuous quantum measurement on a part of the environment to extract information about the system. 

^ ' Here we use the formalism of quantum trajectories |a, S H H to derive a set of conditional master equations which 

Cy describe density matrices inferred by several observers simultaneously performing measurements on different parts 

of the environment. The knowledge a given observer has about the system, inferred from his measurement records, 

leads to his single observer density matrix. One can also consider a super-observer with access to all the records of all 

observers. His knowledge about the system is encapsulated in a super-observer density matrix. Both density matrices 

evolve according to conditional stochastic master equations that will be derived in the next Section. 

The paper is organized as follows. In Section II we derive conditional stochastic master equations for the situation 
C^ ' where multiple observers monitoring the environment of a quantum system try to infer its state from their measurement 
data. Section III discusses correlations between different measurement channels. In Section IV we apply the formalism 
to the quantum Brownian motion at zero temperature, a model for which coherent states are perfect pointer states. 
Observers trying to infer the state of the oscillator, initially prepared in a Schrdinger cat state (|z) -I- | — z))/^/2 of 
large amplitude coherent states, eventually fully agree when they measure in a basis correlated to the pointer states. 
In Sections V and VI we consider the model of resonance fluorescence of a two-level atom, for which pointer states 
do not really exist, and even the most predictable states are quite unpredictable. We discuss correlations between 
single observer state assignments for different measurement schemes, such as photo-detection and homodyne detection. 
Finally, Section VII contains our conclusions. 

II. CONDITIONAL MASTER EQUATIONS FOR SEVERAL OBSERVERS 

Imagine a system S coupled to an environment £. The state of the environment is monitored by a set of detectors 
{T^i} [i = 1, . . . , C) in measurement basis that can be different for different detectors. When the results of these 
measurements are ignored (which in technical terms corresponds to tracing over £ as well as over all the records of 
{Vi}), the reduced density matrix of the system Q{t) evolves according to an "unconditional" master equation (UME) 



dg{t) = dtC g{t) , (1) 

within an infinitesimal time step dt. Here £ is a linear super-operator acting on q. This master equation is "un- 
conditional" in the sense that all the information about the records of {Vi] has been ignored. Note that we are 
assuming that this evolution is Markovian, so that the state of the system at time t + dt only depends on its state at 
time t. In what follows we shall make a still stronger assumption that the state of the fragments of the environment 
prior to the interaction with the system is always the same. In particular, the fragments of the environment that 
have interacted with S do not interact with it again. This quantum Markov approximation is accurate in a typical 
scattering situation. Then there is a natural distinction between input and output fields. The input field is assumed 
to be always described by the same density matrix /x. The output field evolves in a way it must to account for its 
outward propagation. It is entangled with the system, and is eventually measured sufficiently far from the system so 
that the detection does not disturb the system. A spontaneous emission from, say, a two-level atom in free space is a 
particular example of the scattering situation where the input density matrix /z is a vacuum state. 

A. Multiple measurement channels 

Suppose that the measurement records are not ignored but, instead, as is typically the case, they are used to extract 
information about S. Let us further assume that £ includes parts {£i} numbered by the index i (i = 1, . . . , C) and 
that such parts are coupled to the detectors {Vi}. These detectors correspond to different measurement channels. We 
model the measurements by detector Vi as a projection of the detector's state in a measurement basis with outcome 
dNi{t). In the example of spontaneous emission the parts {£i} can be chosen as different directions of photon emission 
monitored by different detectors {Vi}, while dNi{t) may be the number of photons collected by the detector i during 
the infinitesimal time interval from t to t + dt. No matter how much the measurement is delayed or how long it 
takes to decohere the state of the detector Vi, the eventual outcome dNi{t) affects the knowledge about the state of 
the system at the time t when the environment part £i got entangled with S. In general the measurement basis can 
be non-orthogonal and/or over-complete, but in the following, for the sake of simplicity, we restrict ourselves to the 
complete orthogonal case. 

We assume that at i = the state of S was a density matrix g{0). As time goes on, at every time step dt the system 
S is getting entangled with a new fragment of the environment £. After interaction with S some of these fragments 
are not detected, while some other fragments, {£i}, are measured by the detectors {Vi}. After time t the initially 
uncorrelated state oi S + £ +V evolves into an entangled state psev{t) (by V we are denoting the set of detectors 
{Vi}). Ignoring the state of £V we get the unconditional density matrix 

g{t) = Trfp ps£v{t), (2) 

which evolves according to the UME, Eq. (^. Let us call Ma{t) = {Ni{t), . . . , Nc{t)) the set of strings of multiple 
channels measurement records. That is, Afa{t) is a particular history of measurement results {dNi{t')} on all channels 
till time t, and the subscript a denotes a particular set of records. Then, the unconditional density matrix g(t) is a 
sum over all possible sets J\fa (t) , namely 

g{t) = J2 Ti-fP ^A^„(t) PS£-D{t). (3) 

Here Vjsf^[t) = '^JVi(t) ® ■ ■ ■®'PNc{t) projects the state of all detectors according to the measurement records Ni{t){i — 
1, . . . , C). Different branches of the wave- function of S£V are labeled by different sets Na{t) of possible measurement 
records. A branch 7VQ(i) has a probability 

p[Nc.{t)] - Trssv Vu^i^t) Ps£v{t) . (4) 

A hypothetical "super-observer", who knows all the measurement records of all detectors, Ma{t) = 
{Ni{t), . . . , Nc{t)), ascribes to 5 a state given by a multiple observers (or super-observer) conditional density matrix 

p[JVa{t)\ = - — :^ 7- = rrrhvi ' (^ 

Tts Tt£t> V^fU) PS£v{t) PWa{t)\ 



normalized so that Trg p[Afa{t)] = 1- From the point of view of the super-observer the set of records Afa{t) actuaUy 
happens. His description is necessarily probabilistic: at i = he could only calculate probabilities of different outcomes 
(Eq. Q), but could not predict his actual set of outcomes J^a{t)- 

We assign to each measurement channel i an observer i who knows only his own record Ni{t). He ascribes to 5 a 
state given by a single observer conditional density matrix 

pwm - /-y^y)^^^^w (6) 

ir^ ilST, VN,(t) PS£V(t) 

conditioned only on his own record Ni{t). It is easy to check that the single observer density matrix is an average 
over the super-observer density matrices with all records that are not known to the observer i. Indeed, denoting 

Afo '^ ' — {Ni{t), . . . ,Ni(t), . . . , Nc(t)) a set of strings of multiple channel records that contain the particular record 
Ni{t) in channel i, we have 

The probability distribution for the measurement record Ni{t) is given by 

p[N,{t)] = Tvssv VN,it) Ps£v{t) = Y. PiK^^'^l^ (8) 

which is analogous to Eq. Q. 

The issue of compatibility of density matrices ascribed to a system by different observers was first considered by 
Peierls [TJ. He noted that the state assignments of various observers cannot be arbitrarily different, and proposed 
that, in order to avoid contradiction between different state assignments, the product of the corresponding density 
matrices should be non-zero. A second condition put forward in jl2|. namely that the different density matrices should 
commute, was later shown to be too restrictive 13]. Necessary and sufficient condition for compatibility of several 
density matrices turns out to be simple [ij| : Their supports must share at least one state (the support of a density 
matrix is the subspace spanned by all its eigenvectors with non zero eigenvalues). 

In our setting the issue of compatibility of various single observer density matrices p[Ni{t)] is settled very naturally. 
Indeed, the probability p[J^a{t)] defined in Eq. |0J provides a measure of compatibility of different sets of outcomes: 
p[Ma{t)] = when records from different channels in the set Na{t) are mutually contradictory. It is easy to see 
that this condition is equivalent to the pre-requisite compatibility, i.e. the overlap of support [i3j ES 13 • Moreover, 
p[N'a{t)\ quantifies this compatibility, at least in the multiple observer setting we are about to investigate in more 
detail. 

B. Conditional master equations 

The evolution of the unconditional density matrix Q{t) is determined by the unconditional master equation, Eq. 
(^. We now derive master equations that describe the evolution of the super-observer density matrix p[I\fa{t)], and of 
the density matrix p[Ni{t)] of observer i, conditioned on their respective measurement results. Within the Markovian 
approximation, the master equation for the single observer conditional density matrix has the form [a,l3 

ctpWo^it)] = dtL pWc^it)] + MdM^(t) O pWcit)] ■ (9) 

In our case of the super-observer that "single observer" has access to all measurement records J\fa{t). The super- 
operator M.dj\f^(t) conditions p[Na{i)\ on present measurement results dNa[t) — {dNi{t), . . . ,dNc{t)) that are ob- 
tained at time t. It can be written as a sum of super-operators that depend on the measurement results on individual 
channels, ^Ad^fa{t) — MdNi(t) + ■ ■ ■ + M.dNc{.t)- Each super-operator M.dNi(t) takes density operators to density oper- 
ators and depends on the particular measurement strategy implemented in the measurement channel i. Examples of 
measurement strategies are point processes, such as photo-counting of optical fields, and diffusive processes, such as 
homodyne or heterodyne detection of optical fields. Also, MdAfa,{t)® p\^a{t)] is nonlinear in plMait)] and linear in the 



set of measurement results dAfa{t). The nonlinearity comes into play in Eq. ^ when we normalize the density ma- 
trix. The action of this super-operator on p[Afai{t)] is a generalization of the apparent "collapse of the wave- function" 
experienced by the super-observer confined to the branch Afa{t). For a derivation of the super-operators MdNi{t} in 
terms of the projectors V^i^t) see the formalism of operations and effects described in ISlllq. 

When measurement results dJ\fa{t) are ignored, g{t) should follow the UME (Eq. |^). In other words, the sum of 
the super-operator Mdj\fc,it) over all possible strings of measurement results dMait) should vanish, that is: 



Y. Md^^it) Q pWo.it)] = 0. (10) 

dA^a (t) 

We can also write down the master equation for the density matrix p[Ni(t)] that observer z, who knows only his 
own records Ni{t), ascribes to the system S. Using Eas. H7l9|) we obtain 

dp[N,it)] =. dtC p[N,it)] + MdN^^t) p[N^{t)]. (11) 

Indeed, this equation has the same form as that of the super-observer Eq. 0. The super-operator A4dNi(t) depends 
only on the measurement result dNi{t) of observer i, and it is defined as an average over all records of all other 
observers unknown to him. 



MdN.it) &Pmt)]^ ^ ^ ■M^^...(.)0p[<'(*^], (12) 

where dj\fa is any string of multiple channel measurement record that contains the particular measurement results 

dNi{t) on channel i. 

III. CORRELATIONS BETWEEN DIFFERENT MEASUREMENT CHANNELS 

We study correlations between measurement records on different measurement channels, say channels i and j. It is 
clear that in order to do so it is necessary to compare corresponding records Ni(t) and Nj{t). Therefore, we imagine 
there is someone who has access to the records on both channels, or the two observers with access to channels i and 
j communicate with each other and share their measurement records. Whatever the case is, we can think that there 
is a (super)observer who has access to the two measurement channels and whose string of records is {Ni{t), Nj{t)). 
To follow the line of thought of previous sections, we will instead consider the super-observer who has access to all 
measurement channels, and whose string of records contain the particular record Ni{t) on channel i, and the particular 

record Nj{t) on channel j, i.e., his string of records is Afj ^^ ' = (Ni(t), . . . , Ni{t), . . . , Nj{t), . . . , Nc{t)). Here 

the subscript 7 denotes a particular set of multiple channel records that contains records Ni{t) and Nj{t) in channels 
i and j, respectively. 

We define the average relative purity between the states ascribed to S by two such observers i and j as 

Ovit) = E E p[<''*^'^^-'*^]Tr5p[iV.(i)] P[NM 

Ni{t),Nj(t) j^j-i^iW.Njit) 

= Yl p[N^{t),N,{t)]Trsp[N,{t)]p[N,{t)], (13) 

N,{t),N,it) 

where p[Ni{t), Nj{t)] is the joint probability distribution for records Ni{t) and Nj{t), given by 

p[N,{t),N,{t)]^ Y MAC'^*^'"^^'^*^]- (14) 

Ni(t),Nj{t} 



We also introduce the average relative purity Oi (t) between the states of S ascribed by the observer i and the super- 



observer, whose respective measurement records are Ni{t) and Afg , i.e., the particular record Ni{t) on channel i 



is contained in the super-observer string of multiple channel records, Mo '^ ' = {Ni{t), . . . , Ni{t), . . . , Nc{t)). Here /3 
denotes a particular set of super-observer records that contains the record Ni{t) in channel i. We define 

0^it) =T. T. p[M^'^'^]TTsp[Nm pW^'^% (15) 

Using Eqs.|7||Hl) it is easy to check that this average relative purity equals the average purity of the state ascribed to 
S by observer i, and that it is also equal to the autocorrelation Ouit) introduced in Eq. (|13|l . namely 

0,(t) = 0^,(t) = Y. p[N^{t)XYvsp'[N,(t)]. (16) 

In other words, an observer i has no more information about the super-observer's records than the information already 
contained in his own records. 

A better measure of correlations between density matrices is fidelity [13, which is defined as F{pi,p,j) — 
{Tr[(y^Pjy^)^/^]}^. This can be easily calculated in two dimensions: F[pi,pj) = Ti{piPj) + 2(detpidetpj)^/^. 
Unfortunately it is more difhcult to compute for more general cases, and for this reason we will use in the following 
relative purity as a measure of correlations. 

In order to make a quantative study of correlations between different measurement channels it will be convenient 
to work in the framework of stochastic differential master equations. We first simplify our notation: we will denote 
the density matrix ascribed to the system by the super-observer, who has measurement records J\fg , as 

p[AAf (*)] ^ Pit), (17) 

and the density matrix of observer i, whose measurement record is Ni{t), as 

p[N,{t)] = p,{t). (18) 

Also, we will denote the measurement records Af^ '^ ^ of the super-observer, and his measurement results ci/V^ , 
that respectively contain the measurement record Ni{t) and the measurement result dNi{t) on channel i, as 

dM^/'^'^ = dN{t). (19) 

We emphasize again that, in the study of correlations, we consider the situation when the single channel measure- 
ment record Ni{t) is contained in the super-observer records Mp = Mg . Similarly, when we study correlations 
between measurement records on two different channels i and j, we consider the situation when the two single channel 
measurement records Ni{t) and Nj{t) are both contained in the super-observer records M(i = Mo . 

The conditional master equation for the super-observer, that conditions his p{t) on current measurement results 
dM'{t), can be upgraded to a stochastic master equation (SME). A SME is a conditional master equation plus a 
probability distribution P{dM'{t)) for the measurement results dM'{t). This probability distribution can be obtained 
from Eq. Q. P[dJV{t)] is a conditional probability to get current measurement results dM{t) provided that the 
measurement records of the super-observer until time t are Af(t) , that is 

prnm^pidMitwit)]^?^^^^^^, (20) 

where we used Bayes rule. Here p[dM{t),JV{t)] is the joint probability of having measurement records JV{t) until time 
t, and of having measurement results dM{t) at time t. In the Markovian approximation P[dJ\f{t)] depends on the 
records M{t) through the conditional density matrix p{t) , 

P[dMit)] ^ P[dMit)\pit)]. (21) 



The dependence of this probabiHty distribution on the super-observer p{t) leads to correlations between different 
measurement channels we are going to explore using a set of SMEs describing the stochastic evolutions for p{t) and 
Pi (i) . This set of stochastic master equations is given by 

dp{t) = dtCQp{t)+Md^(t)Qp{t), (22) 

dpi{t) = dtCQp^{t)+MdNAt)(^Piit), (23) 

P[dN,it),...,dNaitMt)] = Tr,,.^^.,,,^.,)p,,.ft + dt) ^ ^^^^ 

We refer to this set of equations as a multiple channels stochastic master equation (MCSME). 

We will use the above formalism to address questions regarding correlations between measurements on different 
channels: 

• What is the average correlation between the density matrix of a single observer pi and the super-observer's 
density matrix p? We shall quantify this correlation by the average relative purity Oi (t) defined in Eq. H15|l . As 
a short-hand, we will write it as Oi(t) = Ti pi{t)p{t), where the over line means the weighted average defined in 
Eq. H15|l . This relative purity is a measure of how different, on average, is the knowledge of the observer i from 
the knowledge he would have had if he had access to the records of all the other observers. He cannot know more 
about the state p{t) that the super-observer ascribes to the system that he can infer from his own measurement 
record only. The extracted information can be measured by the average purity of the single observer state 
Oii{t) — Trp|(t) (Eq. (fTE|) '). Even if the super-observer density matrix p{t) had higher average purity, the 
average relative purity Oi (t) would be equal to the single observer purity, Oi (i) — On (i) , as derived in Eq. lfTC|l . 
This equality will be illustrated with several examples in Sections V and VI. Oii{t) is maximal for measurements 
in a basis correlated with the pointer states ;11]. 

• What is the average correlation between different single observer density matrices pi and pj ? We shall quantify 
this correlation by the average relative purity Oij{t) defined in Eq. (|13|l . that as a short-hand we write as 
Oij{t) = Tipi{t)pj{t). In other words, how much do different observers agree about the state of the system? 
In Section IV we will show that for an initial Schrodinger cat state made of large amplitude coherent states 
(coherent states are perfect pointer states in the model of zero temperature quantum Brownian motion), and 
for measurements in a basis of the environment correlated with them, observers will, after an initial transient, 
reach full agreement, Oij{t -^ cx)) = 1. Typically, as seen in the examples of Sections V and VI, the agreement 
is not perfect but it gets better when the observers' measurement basis get closer to those environmental states 
correlated to the pointer basis of the system. For resonance fluorescence from a two-level atom subjected to 
direct photo-detection (see Section V) we find an anti-correlation, Oij < 1/2. Each observer learns something 
about the state of the system but their estimates of the state pi{t) are anti-correlated. The two- level atom is 
very far from being classical and, what is more, photo-detection is very far from being a measurement in a basis 
correlated with the most predictable states. 

IV. QUANTUM BROWNIAN MOTION 

A. Correlation of the outcomes for pointer state measurements 

In this section we consider the well-known model of quantum Brownian motion consisting of a harmonic oscillator 
(the system) interacting with a reservoir of harmonic oscillators (phonons) with a position-position coupling. We will 
restrict ourselves to the case of a zero temperature environment. This model represents a damped harmonic oscillator. 
The self-Hamiltonian for the system is H = oja'^a, where w is the frequency of the oscillator and a, a' are bosonic 
annihilation/creation operators. Imagine that a set of observers perform homodyne detection measurements on the 
environment of phonons so that each of them gains information about the state of the system oscillator. Given the 
set of records of all those observers, the MCSME of the system oscillator is [23 



(a + 7)p(at+r) (2^^ 



dp = dt (apat - l«tap - '-pa^a^ + ^ (dN, - dN.{p) ^ ^^^^^^ ^ ^^^^^^ ^ ^^^^ 

dp, = dt [^ap,a^ ~ -a)ap, ~ -p^a^aj + (dN, - dN,{p,)^ (^_^_2_^L^__2_ _ p^J , (26) 



dN^ip) = r],dt[R^ + Re-"^TTpa + i?e+'"^Trpa^ + Tr pa^ a] . (27) 

These equations for the conditional evolution of the density matrices of the system, written in the interaction picture 
representation, are valid in the rotating wave approximation. Here we use Ito version of stochastic calculus. The first 
terms on the RHS of Eqs. (|25I26|I are of Lindblad form and describe damping and decoherence due to spontaneous 
emission of phonons. We have set the damping coefficient to one. The second (stochastic) terms feed back into the 
master equation information about the state of the system gained by observers. The coefficient 7 = i?exp(i</)) is the 
amplitude of the local oscillator in the homodyne detector |l8| . For simplicity, we are asumming here that all observers 
perform the same kind of homodyne detection, so that the amplitudes Ri and phases 4>i are all equal. We will lift 
this restriction in later examples. The number of phonons detected by observer i in an infinitesimal interval from t 
to t + dt is dNi{t) £ {0, 1}, with an average given by Eq. 1)2 7|l and dNidNj — Sij dNi. The efficiencies rji of different 
detectors can be defined as the fractions of phonons monitored by particular detectors. In the phonodetection limit 
(R — 0), the average detection rate Eq. 127|l is proportional to the average occupation number. Whenever a phonon 
is detected {dNi — 1 for any i) the occupation numbers in p are reduced by one. In the homodyne limit (i? 3> 1) the 
detection rates measure the coherent amplitude Tr/9(e+*'^a + e~''^a^) of the state of the system. 

To illustrate how different observers are gaining information about the system and how correlations between different 
measurement channels arise in the process of continuous measurement we consider superpositions of large amplitude 
coherent states. According to the exact solution '2(T\ coherent states |zo), such that a\zo) = zol^o), decay to the ground 
state like \zq e~*) = \z) without producing any entropy. At T = they are the perfect pointer states of the quantum 
Brownian motion model [13 ■ The decay to the ground state takes place on a time scale of the order of the damping 
rate, which we have set to 1. In a subspace spanned by | + z) and \ — z) a general density matrix is 

Pit) ^^^^\ + z){+z\ + ^—^\- ^){-^\ + C{t)\+ z){-z\ + C*it)\- z){+z\ ■ (28) 

Substitution of this density matrix into Eq. it^ . and subsequent left and right projections on | ± z) ^^ give 
stochastic differential equations for A{t) and C(i). C(i) decays to on a decoherence timescale which in our subspace 
of I ± z) — I ± rexp(i'0)) is given by 1/r^. For initial r ^ 1 this decoherence is much faster than damping and 
it takes place much before the states | ± z) decay to the ground state. In the opposite case of r ^ 1 the states 
I ± z) decay to the ground state before they can be distinguished by the environment. Both limits were considered 
in Ref. |23|. Coherent states with coherent amplitudes ±z that differ less than 1 cannot be distinguished. Here we 
concentrate on the distiguishable case of r 3> 1. In this limit we can self-consistently ignore damping and focus on the 
decoherence and measurement process. In the homodyne limit (i? 3> r), where detection rates are fast as compared to 
the spontaneous emission (decoherence) time, at any given time the correlators for the increments drii = dNi ~ dNi (p) 
are 

drii — , 



dn^duj « 6rj dNi ~ <5y rj^R^dt + 0{R) . (29) 

In this limit the increments can be approximated by drii — ^JrfiR dWi, where dWi^s are gaussian Wiener increments. 



such that dWi = and dWtdWj = Sijdt ja,L3l- After introducing a variable i? as ^4 = tanhB, and translating to 
Stratonovich convention, the super-observer's equation for B reads 

d Fi 

rjta,nhB + J2VniSi{T), (30) 



dr 

i 

where we have defined a new time variable r — Atr^ cos^((/) — iji) and 77 = ^^ rji. Here 9i are stochastic continuous 
functions of time, defined as dWj = Oidt. These stochastic functions are white noises with correlators 



ei{Tl)ej{T2)^5,j5{Ti~T2). (31) 

According to Eq. (jSOJ, B initially performs a random walk driven by the noises but once it diffuses into a positive 
(tanhS = +1) or negative (tanhS = —1) domain, the deterministic force TytanhS takes over and inevitably drives 
B towards positive or negative infinity. After the transient time r ^ I/77, A settles down at A — ±1 which corresponds 
to the pure state | ± z) (see the super-observer's trajectory A{t) in Fig.l). By this time an observer who knows all 




FIG. 1: A single realization of the stochastic trajectories A{t) (thick line), Ai{t) and ^2(1") (thin lines) for rji — 0.7 and 
772 ~ 0.3. The super-observer's A{t) settles at +1 around r « 1, it is followed by ^i(t) after a r-delay ~ 1. ^2(1") after a long 
period of indecision settles down at +1 at r « 5. 

dNiit) can tell whether the system oscillator is in the state \z) or | — z), and attributes to the system the appropriate 
pure state. This happens also when the total efficiency rj is less than one. 

An observer i ascribes to the system a state pi conditioned on his own records dNi (t) only. Since we want to study 
correlations between the measurement records, the evolution of pi is given by the MCSME (Eq. (|26|l ). Taking the 
homodyne R^ r limit in the single observer case we get 



dN, - dN,(p,) = {dN, - dN,[p)^ + {dN,{p) - dN,{pS) "^^^ ^/mR dW, + 2r],rR {A - A,) cosi^Jj ■ 



(32) 



Just as for the case of the super-observer (Eq. H25|) ). substitution of the ansatz Eq. H28I) into Eq. H26I) and neglecting 
any 0{1/R) terms, leads to the equation for the single observer Bi, 



= T]i tanh Bi 

dr 



[rii tanh B — iji tanh Bi + ^/rfi 6i] = i]i tanh B 



(33) 



where Ai — tanhB^. The terms in the square brackets come from the stochastic term in Eq. 1)26(1 . Note that the 
super-observer's B appears in the evolution equation of the Bi associated to the single measurement channel i. This 
reflects the fact that the single channel and multiple channels measurement results are correlated in the MCSME (Eq. 

Let us now study how the evolutions of A according to the super-observer and single observers are correlated. On 
the one hand, according to Eq. H3U|) . the super-observer evolution settles A = tanhi? at ±1 after the transient time 
r ~ 1/77. On the other hand, the single observer evolution is given by Eq. (|33|) . and correlations between the two 
evolutions enter through the first term in the most right hand side of that equation containing the super-observer 
A = tanhS. Once A = tanhS = ±1 is chosen, the deterministic drift term rji tanh B on the RHS of Eq. 1(33(1 
will inevitably force Ai = tanhi?i to make the same choice after the longer transient time t ^ l/j/j. Eventually all 
observers will settle down a.t A — Ai — ±1, and the average relative purities will be equal to one, Oi = Oij = 1 (see 
the single realizations for two measurement channels z = 1, 2 in Fig.l). 

In our example the observers finally find out which of the two coherent states is the state of the system. It is 
possible because the initial coherent states have large amplitudes ±z with \z\ ^ 1 so that the decoherence time is 
much shorter than the spontaneous emission time. In the opposite regime of |z| <C 1 the decoherence time is longer 
than the spontaneous emission time, and the observers will not find out the state before it decays to the vacuum state. 



B. Independence of the outcome distributions 



To have a better feeling of the multiple channels stochastic master equation (MCSME) formalism we consider the 
following example. There are two observers z = 1,2. The measurements by observer 2 affect the evolution of the 
super-observer's density matrix p. Since the environment is also monitored by observer 1, in principle observer 1 



may be able to identify perturbations of p produced by measurements of observer 2 and realize that there is another 
observer monitoring the system. In Section II we gave a general argument that, as a direct consequence of the quantum 
Markov approximation, observer 2 cannot find out if there is another observer. Here we present a simple calculation 
which illustrates this fact in our example of superposition of coherent states. 

To begin with, note that Eq. (I30|) is equivalent to the following Fokker-Planck equation for the probability distri- 
bution P(t, B) for B at time r 23], 



— ^ = -T^tanh(B)P+ - — 
J] dr dB ^ ' 2 dB 



r— = -7n?tanh(B)P+-,^P . (34) 



We can compare the following two situations: 

1) Observer 1 is the only observer or 772 = 0. His probability distribution evolves according to Eq. (|34() . The initial 
condition P(0, Si) = 5{Bi) leads to the solution 

P'"(r.BO^ '^^' K""^^ K (35) 

V27r77ir ^/27rriiT 

2) There is an observer 2 with 772 3> ryi- In this limit, where the perturbations by observer 2 are the strongest, 
one is most likely to suspect that the less efficient observer 1 could find out about the more efficient observer 2. The 
evolution of B{t) is mainly conditioned upon the measurements of observer 2. The multiple observer A settles at ±1 
on a time-scale 1/r/ which is much faster than the time l/?yi observer 1 needs to find out about the system. The state 
of the system is settled without any influence of measurements of observer 1. Suppose that, with the probability 1/2, 
the multiple observer state tanhS = +1 was chosen. Bi evolves according to Eq. H33(l with a fixed tanhS = +1. 
The probability distribution for Bi is 

P+(r,Bi) = ^^^ ^ . (36) 

y/2'KrilT 

Also with the probability 1/2, tanhP — —1 can be chosen. Now tanhP ~ —1 is fixed in Eq. (|33l) and the probability 
distribution is 

P-{t,B,) = ^ ^^ ''-' ^ . (37) 

As we do not know which super-observer's state will be chosen, the two probabilities add to give P'^^^ — P^ + P_. It 
is easy to check that P'^^^ — P^^-' in Eq. (|35|) . The probability distributions coincide, so observer 1 cannot find out 
if there is any observer 2 even if observer 1 detects just 1% of phonons and the other more efficient observer detects 
99% or almost all phonons. 

V. TWO-LEVEL ATOM: DIRECT PHOTO-DETECTION 

We want to contrast the quantum Brownian motion model with an example of a system with a small Hilbert space, 
such as a driven two-level atom coupled to the radiation field, for which we do not expect perfect pointer states. In 
Appendix A we derive the MCSME for a two-level atom driven by a laser beam with frequency to and whose emitted 
radiation is subjected to photo-detection. It takes the form 

dp = -^ dt [cua.,,p] + dt (cpc^ - ^c^cp - ^pc^c) + ^ {^dN, - dN^ip)) ( rj.T„ - p) , (38) 

dp, = -i dt [iuax, Pi] + dt i cp,c) - -c^cpi - 2'°*'^^'^) '^ (^^* ~ dNi{pi)j i ' ^, - Pij , (39) 

dMjjJ) = r],dtTr[pc^ c] . (40) 
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The density matrix p of the atom is a 2 x 2 matrix 

P = -[I + xa^ + yo-y + za^] . (41) 

The lowering operator is c = {ax — *f]^)/2, and the number of photons detected in channel i between t and t + dt is 
dNi e {0, 1} with an average proportional to the occupation number of the atom, see Eq. (gU)), and dNidNj — SijdNi. 
Following each detection of a photon (any dN^ = 1), the atom is known to be in the ground state (the —1 eigenstate 
of ffz), from where it is excited again by a laser beam through the Hamiltonian term coax- The efficiency rji of the 
detector used by observer i is the fraction of photons which are detected by him. 

When u! ^ 1 the most predictable states of the two-level atom are a^ eigenstates, i.e., they are determined by 
the Hamiltonian cuax describing the excitations via the laser beam |2J|. These states are far from perfect since they 
have a nonzero initial rate of purity loss. Moreover, while eigenstates of (Jx are most predictable, they are not the 
most effective in making an imprint on the environment (as real pointer states should be 4. 9. lOj. In particular, the 
environment-system Hamiltonian does not preserve them. As a consequence, we do not expect agreement between 
observers even if they are measuring in a basis of the environment correlated to the cra;-eigenstates of the atom. Direct 
photo-detection is a way to find out if the atom is in the ground state. This state is complementary to the most 
predictable states. That is why we expect the relative purity between observers to be very poor. In fact we will find 
any two observers to be anti-correlated, Oij < 1/2. 

A. The w > 1 limit 



For w ^ 1 Eas. (|38l39|l can be solved rigorously. Suppose that no photons are detected for a certain period of time, 
dNi{t) = 0. During this time the density matrix p in Eq. (|41|l evolves according to the deterministic part of Eq. (J3SJ). 
The unitary self-evolution with the Hamiltonian luGx is mixing y and z with the frequency 2lj. It is convenient to use 
the interaction picture, where 

y = 2/int COS 2cjt - Zint sin 2cji, (42) 

Z — 2/int sin 2Ljt + Zint COS 2ujt, 

and the variation in time of a;int, J/int, -^int is slow as compared to lj. When we substitute the density matrix Eq. 141|) 
into the deterministic part of Eq. H38() . use the interaction picture, and average over one period of oscillation with 
frequency lo, we will obtain the following equations |29| 

dzint 1 



dt 2 

dy\nt _ _3 

di ~ 4 

dz\^t ^ 3 

dt ~ ^4 



(l-'7)2/int, (43) 

(1 -77)Zi„t. 



Every time a photon is detected the super-observer density matrix p is projected to the ground state. All the 
information about the previous evolution of p{t) is forgotten. Suppose that a detection took place at time i = 0. 
Just after the detection the initial conditions are a;(0+) = 0,j/(0+) = and z(0+) = —1. Before the next detection 
happens, x,y,z evolve according to Eas. H42l43|l . 



(44) 



where t is the time elapsed since the last photo-detection. This solution is valid until the next detection takes place. 
The next detection at t = td will bring p to the ground state again, from where the system will be excited according 



mm^ 


0, 








X{t) = 




Y{t) = 


3 

e * 


(1^ 


-v)t 


sin 2wi, 


Z{t) = 


— e^ 


-K 


'J-1 


)*cos2a 



11 

to a; = X{t — td), y = Y{t — td), z = Z{t — td)- The probability that an observer i will detect a photon between t and 
t + dt after the last detection by any observer is 

1 + Z{t) 



dN,{p{t)) = dtj^, ^. (45) 

The above argument can also be applied to Eq. H39|l . Every time an observer i detects a photon his state pi jumps 
to the ground state, from where it is excited according to 

X,it) = 0, 

Yi{t) = e-t(i-''')*sin2wt, (46) 

Z,{t) = -e-3(i-').)tcos2wi. 

The time t here is the time since the last detection by the observer i. 

B. Distribution of waiting times 

In this example we shall see again that a single observer cannot find out if there is any other observer. We will 
consider just two observers i = 1,2 and we will derive the distribution of waiting times (times between subsequent 
detections) for observer 1. We will show that this distribution does not depend on 772 so it is not sensitive to the 
presence or absence of any observer 2. Any higher order correlations between detection times can be expressed by 
this distribution of waiting times because every time a photon is detected by observer 1, the atomic state goes down 
to the ground state so that any history before the detection does not affect evolution that follows that detection. The 
distribution of waiting times contains all the information observer 1 can possibly extract from his measurements. 

Suppose that observer 1 detects a photon at time i = 0. What is the probability wi(t) that he will detect the next 
photon at time t = t 1 If observer 1 were the only observer, so that 771 = rj, then the answer would be 

wi[T)^\rii ^ je Jo Ri _e 2 . (47) 

The first factor is the average detection rate Eq. (|45|l . and the second one is the probability that no photon is detected 
between and t. As it should be, wi(t) is normalized to unity. To obtain the final expression for wi{t) in Eq. H47|) 
we have neglected all terms which vanish for a; ^ 1 as well as fast oscillating terms ^ cos 2uot. 

If there is a second observer, then the detection rate of observer i depends not on Za{t) but on z{t) (see Eq. (03)- 
In general there may be n — 0, 1, . . . ,cx) detections by observer 2 between and r. Every time there is a detection 
by observer 2 at i = ij, [j = 1, . . . ,n), z{t) jumps down to —1. For tj < t < ij+i it evolves as z{t) = Z{t — tj). 
The probability that there is no detection by observer 1 between times and r, given that there are n detections by 
observer 2 at the times ti, . . . ,t„, is given by 

D„(ii, . . . , t„,T) = q2{ti)q2{t2 -h)... q2{tn - i„_i)e-t(--*")^ (48) 

where 92 (t) is distribution of waiting times for observer 2 given that there are no detections by 1, 

(?2(t)"«' 6-^^^(1-6-^(1-")^ cos 2C.T), (49) 

and the factor e-a^*"*") is the probability that no detections by any observer take place between i„ and t. The 
distribution of waiting times for observer 1, averaged over detections by observer 2, is given by Z?„ multiplied by the 
detection rate of observer 1 at r, and averaged over all possible n and ti, . . . , t„. Therefore the final expression for 
the waiting time distribution /wait (t) for observer 1 in the presence of detections by observer 2 is given by 



O -,-7- -,-;- -,-7- 

^ / dti dt2... dUDn{tu. 

^qJo Jti Jt„-i 



U.itiT) = J2 dh dt2... di„i?„(il,...,i„,r)|-[l + ^(T-t„)]"«'|-6-'^^ (50) 
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where, again, we have neglected terms which vanish for cj ^ 1 and any fast oscillating terms. In Appendix B we show 
how to obtain this last formula. We conclude that the distribution of waiting times for observer 1 in the presence 
of detections by observer 2 (Eq. (jSOJ) is the same as that for no observer 2 present (Eq. (@Z|)). The distribution of 
waiting times for observer 1 is not sensitive to observer 2. 

C. Average relative purity betAveen pi and p 

Let us now find out how much does a given observer, say i — \, know about the state of the super-observer. To this 
end we will calculate the average relative purity between the single observer and the super-observer density matrices, 
Oi = Trpip. Imagine the following situation. Take an arbitrary instant of time t — and call o" {t — 0) — Trpi (r — 
0)p{t = 0) the relative purity given that the last detection of observer 1 took place at r = — i and there were n 
detections by observer 2 between t — —t and r = 0. Then the average relative purity Oi evaluated at the time 
r = will be equal to the <- average (i.e., average over all possible initial times of detection by 1) of the relative 
purity o]^ {t = 0) given that there were no detections by observer 1 between t — —t and r = and averaged over 
all the possible numbers n of detections by observer 2 and his detection times ii, . . . ,t„. For the sake of clarity, we 
now shift the time origin as r ^ r -I- 1, so that the last detection of 1 took place at time and we are interested in 
evaluating Oi at time t. The unnormalized probability distribution for no detections by observer 1 between and i, 
and n detections by observer 2 at the times ii, . . . , t„ is D„(ii, . . . ,tn,t), given by Eq. 148() . The normalizing factor 
for this distribution is 

"1 = 1]/ dt dh dt2... di„D„(ii,...,t„,t)"« — . (51) 

Given that the last detection by observer 1 took place at time and the last detection by any observer happened at 
time i„, the relevant relative purity is 

o["\t) = Trpi{t)p{t) = i + L-i(^-m)t^-U^-v){t-U) cos2u;i„ . (52) 

This relative purity, when averaged with the probabihty distribution H48|) . gives 



o<3 „oo pt pt 

Oi = "-r^ X! / '^w ^^1 / '^^^ 

„=n "'o ^0 Jti 



dtn Dn{tl, ...,tn,t) 0^r\t) =2 + 2(3 ^ V) ' ^^^^ 



where we have neglected all 0{1/uj) terms. 

On the other hand, the average purity gained by observer 1 can be calculated as follows. According to Eas. (|46|l . 
the purity at the time t after the last detection is oii(t) — Trp^(i) — -^ + ^exp[— 1(1 — rii)t]. The probability that 
there was no detection between and t is exp(— ryit/2). Therefore, the average purity On is oii(i) averaged over t: 

_ /o°°rf^e-^*oii(t) _ 1 r,, 

^"" Tdie-^* -2 + 2(3-2,1)- ^^'^ 

As expected from Eq. (|16|l . On coincides with Oi (see Fig. 2). Let us now comment on the limiting cases rji = and 
r]i = 1. In the former case we get On = 0.5, that corresponds to no information gain by the observer (pi is maximally 
mixed). In the latter case we get On = 1, that is maximal gain of information, and pi is pure. 

D. Average relative purity Oij 



The average relative purity O12 = Trpip2 has contributions from the following two situations: 

(1) The last detection by observer 1 took place at time 0. Between times and t there were n > 1 detections 

by observer 2 at the times ti, . . . ,t„. The last detection before t was made by observer 2 at time t„. According to 

Eqs.|03, the relative purity at t is 
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FIG. 2: Average purity On — 1/2 and average relative purity Oi — 1/2 for an observer performing photo-detection measurements. 
The initial condition is maximal lack of knowledge, i.e. p(i = 0) = pi(i = 0) = 7/2. The super-observer's efficiency is 77 = 0.6 
and the single observer one is 771 = 0.5. According to Eas. 1531 1541 the asymptotic value is Oi(oo) = 0\\{oq) — 0.125. The 
stochastic trajectories are an average over 256 single realizations. The small discrepancy between Oi and On in the figure is 
an artifact of the finite number of realizations used for calculating the averages. 



o('2')(t) = Trpi(i)p2(t) 



i + ie-^(i-''i)*e-i(i-''^)(*-*") cos2cc>i„ 



(55) 



The normalizing factor for the probability distribution is 



ni2 = 2J / '^M ^^1 / ^^2 

n=l' " 











tl 



t^>l 



2V2 



dtri J-Jji 1^1 , . . . , in-, t j - , , , 



(56) 



and the averaged relative purity is 



O12 



H2YI 



dt I dti / dt2 
Ja Jti 



dtnDn{ti,...,tn,t)o^T^\t) = -- 



?yi(»7i +m) 



772(6 - 771 - r72)(7 - 772 - 4?7i) ■ 



(57) 



(2) The last detection before t was made by observer 1 instead of observer 2, as in the case (1). The description 
of this second situation is the same as above, except that observers 1 and 2 are interchanged. In particular, the final 
result for the relative purity reads 



?72(?yi +V2) 



_(2) t 
^^ 2 771(6- ?7i - 772) (7 - 771 - 4772) 



(58) 



In general 771 ^ rj2 and the two situations are not equally likely. Let us call p'^^ the probability that case (1) happens; 
clearly for case (2) we have p^^^ = 1— p^^^. The probability p*^^^ is given by p^^-* =7112/^1 = V2/{vi+V2)- The relative 
purity averaged over the two situations is then 



O12 = 



V2 



.(1) 



m 



012 



m + m m + m 

1 771772 [6 - 2(771 +772)] 

2 2(6 - 771 -772)(3-277i)(3- 2772)' 



(59) 



In figures 3 and 4 we show simulations of the time evolution of the relative purity O12 for the case 771 = 772 and t/i y^ 772. 

Note that the average relative purity is manifestly less than 1/2: the single observer states pi and p2 are anti- 

correlated. The reason for this anti-correlation can be explained as follows. Suppose that the states pi-,P2,P are 

initially fully correlated (i.e., relative purity equal to one). Observer 1 is most likely to have a detection when the 
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FIG. 3: Average relative purity O12 — 1/2 between two observers performing photo-detection measurements. Their initial 
condition is maximal lack of knowledge, i.e. pi(i = 0) — pi{t = 0) = 1/2. The efficiencies are ryi = 772 = 0.5. According to Eq. 
1591 1 the asymptotic value of the relative purity is -0.025. The stochastic trajectory is an average over 256 single realizations. 
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time (t) 

FIG. 4: Average relative purity O12 -~ 1/2 between two observers performing photo-detection measurements. Their initial 
condition is maximal lack of knowledge, i.e. pi{t — 0) — pi{t = 0) = 1/2. The efficiencies are -qi = 0.7 and 772 ~ 0.3. According 
to Eq. 15911 the asymptotic value of the relative purity is -0.022. The stochastic trajectory is an average over 256 single 
realizations. 

super-observer's state is excited (z w +1). The hypothetical positive correlation means that when z w -1-1, then 
also zi ~ +1 and Z2 ~ +1. Suppose that a detection by observer 1 happens. The super-observer z and the single 
observer zi jump down to —1. The observer 2 has no clue that there was a detection by observer 1. What is more, 
the super-observer z is close to —1 so observer 2 cannot detect a photon and jump to Z2 = — 1. His Z2 remains close 
to -|-1. Just after the detection the product z zi > but the product zi Z2 < 0. This mechanism cannot make 
Oi < 1/2 but it can and it does make O12 < 1/2. 

We have solved exactly the problem of correlations between multiple measurement channels in the limit of w 3> 1 . 
This limit is sufficient to illustrate our ideas. However, the exact solution for arbitrary u of the resonance fluorescence 
problem in Ref. 26] suggests that, with some extra work, our formulas for average relative purities can be generalized 
exactly to arbitrary lo. 



VI. TWO-LEVEL ATOM: HOMODYNE DETECTION 



As we saw in the previous section, direct photo-detection is a way to find out if the atom is in the ground state. 
One can also measure different quadratures of the two-level atom by performing homodyne detection on the radiation 
emitted from it |l3 ■ In general, it is possible to measure the expectation value of the operator {x cos 4> — y sin (f) , 
where is the phase of the local oscillator in the homodyne detector. This kind of measurement tends to localize the 
state of the atom around the eigenstates of the operator [ax cos(p — ay sin0). The MCSME is (see Appendix A) 



dp 



'i dt [ojax, p] + dt ( cpc^ — —c^cp — T:pc c ) + 



J2 {dN, - dN,{p) 



(c + 7,)p(ct + 7*) 



15 



(60) 



dpi = —i dt [ujax, Pi] + dt cp. 



Tr[(c + 7.)p(ct+7*)] 



dNiip) = r]idt[RJ + i?,e-"^-Trpc + i?ie+"^'Trpct + Txpch 



(62) 



Here 7.^ = Ri exp(i0,j) is the complex amplitude of the local oscillator of the detector i. We will eventually take the limit 
Ri — > 00. We allow each observer to have his own homodyne phase 0,;, so that they can measure different quadratures, 
i.e., they measure non-commuting observables (a related experimental realization of measurements of non-commuting 
observables in two channels in cavity QED was carried out in [23). The detector currents are proportional to Eq. 
(|62|l . The case 0i = corresponds to measurement of the x— quadrature and (fn = 7r/2 to y— quadrature. The large 
Ri limit of Eq. (EOJ is 



dp 



-i dt [uj<Jx, p] + dt I cpc 



-c'cp 



■Y,VvldW,[cpe-"^^+pc^e 



U+i<l>i 



-pc^c 



pi:ricpe-">'^ +pc^e+"*'^) 



(63) 



where dWi are gaussian Wiener increments such that dWi — and dWidWj = Sijdt. To derive the large Ri limit of 

Eq. (inH) we first split dNi - dN~{pi) = (dN.^ - dN~{p)] + CdMjjjj ~ dN.i{pi)\ . In the large Ri limit the first term 

is proportional to Riy/rfidWi, and the second term is proportional to RiriidtTi[c{p — pi)e^^'^^ + (p — pi)c'e'^^'^']. The 
large Ri limit of Eq. l|HT)) reads 



dpi = -i dt {uJUx-, Pi] + dt ( cpiC^ - -c'^cpi - -piC'c 



(64) 



+ [V^ dW, + T^dt Tr {c{p - p,)e"'*- + [p - p.)c^e+'^')] [(cp,:e"'*' + P^c^e+'^') - pTr (cp^e"^"*' + ^0^6+^^')] 



A. Average relative purity between pi and p 



Unfortunately it is not 
For small values of these 
in powers of 77^. For 77^ = 
has a stationary solution 
matrix is perturbed from 
Tji. We expand p — pss + 

We expand a; as a; = x^^' 

used for y and z. To first order in ry^ Eq. H63|l reads 



possible to find analytic solutions to the above equations for all values of the efficiencies T]i . 
efficiencies it is possible to work out various relative purities by a perturbative expansion 
0, the conditional master equation (Eq. (|63|l ) is the unconditional master equation, which 
Pss- In the limit w ^ 1 it is equal to pss — -f/2 or Xss = Vss = Zss = 0. The full density 
this stationary state by the noises dWi , and the magnitude of the perturbation grows with 
(5p, the last term containing those perturbations. Let us write 5p — {xux + yjy -f- zaz)/'2-- 



.(2) 



where x^^' is of order 



Vt 



1/2 



^"^1 is of order ry^ , etc. Similar expansions are 



These equations have a solution 



dt 

dy^ 
dt 

dz^ 
dt 



X' 

2 



y^ V^^i C0S(/)i, 



2^ 



--y^') -2ujz'^^^ -^y/Tf^di sine 



z(^) + 2Luy 



(1) 



(65) 



xW(t) 
y(i)(t)±iz(i)(t) 



i 



±^z^\i) 



(66) 
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FIG. 5: Average purity On — 1/2 and average relative purity Oi — 1/2 for an observer performing homodyne measurements. 
The efficiency is ryi =0.1 and the homodyne phase is (jf>i = 0. According to Eq. H7()|l . which is vaUd for smaU efficiencies, the 
asymptotic value of the average relative purity and average purity is 0\ — 1/2 = On — 1/2 — 0.05. In the scale of the figure 
Oi and On practically coincide. The stochastic trajectories are an average over 256 single realizations. 



where 



4^\t) - Vm COS0, / dT e-^(*-") ^,(t), 

J — OO 



(67) 



To leading order in rj^ 's the single observer equation (|61|l is 



dx. 



(1) 



dt 



dy. 



(1) 



dt 

d^ 

dt 






-^+^-^- 



These equations are solved by the already introduced x\ ,y 
Oi = Tipip is 



(1) „,(i) Ji) 



(68) 



To leading order in rji's the relative purity 



O, 



1 1 

2 + 2 



5(l)a;(l)+y(l)y(l)+ ^(1)^(1) 



1 1- 

2 + 2 



•^i -^i ^ tli iJi T" -ij ^j 



A straightforward calculation leads to the following stationary average relative purity 



Oj = - + 77j 



- cos 
2 



— sni 
3 



(69) 



(70) 



As we can see from Eq. (I69|l the average relative purity coincides with the average purity On. The latter is the 
highest for measurement basis correlated to the pointer states basis of the system, i.e., when (/>,; — 0. Through this 
measurement one can find out most about the system. In figures 5 and 6 we plot the average relative purity Oi and 
the average purity On for different values of the efficiencies and homodyne phases. 



B. Average relative purity Oi 



The average relative purity O12 — Ti'PiP2 is zero to leading order in r/i. To get a nonzero average relative purity we 
have to go one step further in the perturbative expansion for x, y, and z. The equations for the second order terms 
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FIG. 6: Average purity On — 1/2 and average relative purity Oi — 1/2 for an observer performing homodyne measurements. 
The efficiency is rji = 0.5 and the homodyne phase is (f)i = 0. We do not expect Eq. 17UII to hold for such a big efficiency. The 
stochastic trajectories are an average over 256 single realizations. The small discrepancy between Oi and On in the figure is 
an artifact of the finite number of trajectories used to calculate the averages. 

that follow from the single observer equation Eq. (|64|l are 



dx 



(2) 



1 (2) 



— = - -xl"' + 7]i cos^ (j)ixy' - r]i sin 0i cos (jyiyy^ , 



dyl 



(2) 



dt 



,(2) 



.(2) 



.2w..„(l) 



dz. 



(2) 



dt 



-yl ' -2ujzl ' + ?]i sin 0iy]- - rji sin^i cos^iX^ 



(1) 



Formal solutions of these equations are 



(71) 



x?\t) = V^ dT e-^(*-^) (cos^ <P,xf\T) - sincj), cos(/),yf V)) , 

J —OO 

yf\t) ± izf )(t) = r,J dT e-|(*--)±2-(*-) (sin^ ^.yf V) - s^'/'.: cos</),xf V)) • 
To the first non-vanishing order in r]i the average relative purity is 



(72) 



1 1. 



Oi2 = I + ^[a;ia;2 + yiy2 + Z1Z2] = 2 + 2 



^(1)^(2) , ,,(i),,(2) I ^(1)^(2) 



^(2)^(1) , ,,(2) (1) , (2) (1) 
•''1 ■''2 ~r i/l i/2 ^ ^1 ''2 



We evaluate this expression in Appendix C. The result is 



Ol2 = 



?yi//2 



cos^ 01 cos^ '?^2 + 7: sin^ <f)i siv? <j)2 



(73) 



(74) 



The average relative purity is maximized when both observers perform x-measurements (01 =02 = 0). We verified 
this formula by numerical simulations using rji = r]2 = 0.01. Below, in figures 7 and 8, we plot the average relative 
purity O12 for different sets of homodyne phases and efficiencies r/i = 772 = 0-1. These efficiencies are beyond the 
range of validity of Eq. 1741) . 



VII. CONCLUDING REMARKS 



Let us summarize the new results contained in this paper. We have studied continuous quantum measurement 
with several observers and we have demonstrated that it reduces to the "single observer" case. The key problem 
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FIG. 7: Average relative purity O12 — 1/2 between two observers performing homodyne measurements. The efficiencies are 
rji — Tj2 — 0.1 and the homodyne phases are (j)! = 4)2 = 0. The stochastic trajectory is an average over 256 single reahzations. 
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FIG. 8: Average relative purity O12 — 1/2 between two observers performing homodyne measurements. The efficiencies are 
rji — ri2 — 0.1 and the homodyne phases are (j)i = (f>2 ~ 7r/2. The stochastic trajectory is an average over 256 single realizations. 



of consistency of the sets of data acquired by different observers is then reduced to the probability that a given 
combination of data sets will be ever detected by the super-observer. We have applied the formalism to several 
examples of quantum optics as well as to quantum Brownian motion. Observers gain information about the state of 
the system from their measurement records. We have shown that observers gain most information about the system 
and they agree the most when they measure in environmental basis most correlated to the pointer basis of the system. 
Several questions regarding correlations between measurement records of different observers were posed. We have 
shown that the problem of consistency of sets of data acquired by different observers is reduced to the probability that 
a given combination of data sets will ever be detected by the super-observer. We have introduced average relative 
purity to study correlations between measurement records of different observers. For the model of zero temperature 
quantum Brownian motion (which is equivalent to the model of a damped harmonic oscillator) coherent states are 
perfect pointer states: The solution to Eq. H25|) for an initial coherent state remains pure and it is just a coherent 
state with decaying amplitude. We have shown that for an initial Schrodinger cat state {\z) + \ — z))/\/2 made of 
large amplitude coherent states, records of different observers performing measurements on the environment in a basis 
correlated with the pointer basis will eventually fully agree (as shown in Figure 1), and the average relative purity 
will be equal to one. For the case when the most predictable states exist, but are not very predictable and are not 
imprinted on the environment (and, in particular, do not commute with the interaction Hamiltonian), such as the 
model of two-level atom resonance fluorescence, the agreement between observers' guess of the state of the system 
may only be partial, and it is even possible to obtain anti-correlation between measurement records, as in the case of 
photo-detection. 
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APPENDIX A: DERIVATION OF THE TWO-LEVEL ATOM PHOTO- AND HOMODYNE DETECTION 
MASTER EQUATIONS FOR MULTIPLE MEASUREMENT CHANNELS 

Let us assume a two-level atom that interacts with the electromagnetic field, which we shall consider as the 
environment. We will split this environment into different parts i, each of which has associated a detector i. For 
example, i may denote different photon wave vectors. In the rotating wave approximation, the dipole interaction 
between the atom and the electromagnetic field is 

V{t)^tY,{blc-c^h), (Al) 

i 

where bi and c are annihilation operators for photons and the atom, respectively. At every instant of time t, a new 
part of the environment is interacting with the system. Indeed, a localized photon wave packet arrives at the atom, 
interacts with it, and then flies away. Subsequently, a new wave packet performs the same process. Imagine that at a 
given instant of time t the combined state of the atom and the field is R{t) = pit) (E) /i, where p is the density matrix 
for the atom and p is that for the field, which we asumme to be in vacuum, p = ®i ^ |0) (0|i. This series of idealizations 
are called the quantum Markov approximation. 

The evolution operator for a time interval dt is U{t^t + dt) — exp[^,-((i_B|c — c^dBi)], where dBi{t) — bi{t)dt has 
commutation relations 

[dB,{t),dB]{t)]=S,jdt (A2) 

that follow from the (singular) commutation relations [dbi{t), dbUt')] — SijS{t — t'). The above commutation relation 
is of order dt instead of dt'^ [7|, as one might have naively expected. For this reason an expansion to first order in 
dt of the evolution operator requires a second order expansion in terms of dBi and dBj . When ones discards all the 
information contained in the environment (which is then traced out) one gets an unconditional master equation for 
the system 

dp = —idt[u}CFx, p] + dt\^{cpc^ — ■;^c cp— -pc)c). (A3) 

i 

The sum over i just re-scales the spontaneous emission rate of the atom. In the following we shall absorb that rescaling 
in a redefinition of time and set the spontaneous emission rate to one. 

If the measurements on the environment are not ignored but kept, the evolution of the system is conditioned upon 
them. In the case of photo-detection, for most of the time intervals no photons are detected. In this case of null 
results the density matrix of the system evolves according to 

d/Ozcro(0 = O,cro/5(i) - p(t)Tr(acro/5(i)), (A4) 

Ozcro/0(t) = dt 



Uz[i/,p]-i^{ctc,p}j, (A5) 



which is so constructed as to conserve the trace of p under the time evolution. When a photon is measured by any of 
the detectors, the system discontinuously jumps to the ground state of the atom 

OoncPit) = riicpc^dt. (A7) 
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Here the increments dNi € {0, 1} are dichotomic stochastic processes with averages dNi{p) = Ty^diTr [pc^c], -qi denotes 
the fraction of the environment measured by detector i, and fi is such that two conditions must be satisfied: f) 
when dp = dp^cm + dponc is averaged over all records i, it must reduce to the unconditional master equation, and 
;2) Tr[(ipzoro + dponc] — 0. ft then follows that /,; = —riipTT(OzcioP)/dNi. Finally we get the super-observer master 
equation for photo-detection 



dp — —idt [cjax, p] + dt I cpc^ — —c'cp — -;::pc^c j + ^^ ( dNi — dNi[p) 



( cpc? 



\ Tr [cpct 



(A8) 



The super-observer unconditional master equation (jA3|) is invariant under the transformation c ^ c + 7.^ and 
H ^r H — {i/'2)'^i{'y*c — 7iC^), where 7, is a complex number |3|- This symmetry is helpful for deriving other 
unravelings of the unconditional master equation, for example the one corresponding to homodyne detection. In this 
case 7j represents the coherent amplitude of the classical field of the local oscillator i. Introducing this symmetry into 
the photo-detector master equation one immediately obtains the homodyne master equation Eq. (|6U|I . 



APPENDIX B: CALCULATION OF THE DISTRIBUTION OF WAITING TIMES 

In this appendix we calculate the distribution of waiting times /wait (''') for the model of resonance fluorescence from 
a two-level atom subjected to direct photo-detection. It is given by Eq. H50|) . 



OC ^-7- p'j' ^-7- 

/wait(T) = ^ / dti dt2... dtn Dnih,. ..,tn,T)-;^[l + Z{t -tn) 



(Bl) 



where we recall that 



and that I?„ is 



Z{t - i„) = _e-3(i-'')(--*") cos2a;(T - t„), 



A.(ii,...,i„,r)=e-^(-*")nfe-i(*^--^) 



e 4\ 



' cos2u!{tj — tj-i] 



j=i 



(B2) 



(B3) 



where tg = is the time of the last detection by observer 1. Inserting this equation into the previous one, we see that 
when doing the n time integrals only two terms will survive: one that steins from the product of all the I's in Z?„, 
and another coming from the products of all the cosines (which will therefore contain factors of the form cos'^{2ujtj)). 
All other terms in the expansion of the product in _D„ will vanish upon integration. In the lu ^ 1 limit we can replace 
cos'^{2ujtj) by 1/2. We then get 



/wait (t) 






n=0 

00 r 



= T^ ^E 



n=0 






o 



dti . . . dtn 

Jt„-i 

i_(mtV 1 

n! V 2 
1 



-^)" 6-^(1-")" cos(2u;r)' 



_Lf!^l __Lf_!Ml p-f(i-'?)- 



n 



cos{2llit) 



(B4) 



APPENDIX C: CALCULATION OF THE RELATIVE PURITY FOR THE TWO-LEVEL ATOM WITH 

HOMODYNE DETECTION 



In this appendix we derive Eq. H74|l for the stationary value of the average relative purity between two measurement 
channels for the model of resonance fluorescence from a two-level atom subjected to homodyne detection. 
We must calculate the different terms of Eq. (|73|) . Using Eq. (|72ll we have 
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4'\t)^fit) = Vj f rfre-^(*"") (cos^ 0,xf ^(i)x|'^(r) - sin0, cos<^,a;f ^(i)yf ^(r)) . (CI) 



Using Eq. (|^ and that dW-idWj = 5ijdt, it is easy to show that x] {t)yi (r) = 0{1/uj), so wc can discard that 
term in the previous equation. Also, xi {t)xi (r) = rji cos^ (pi exp(— (i — t)/2). Hence 



2^1 {i)Xj {t) = ^j'yj cos^ 4>i cos^ 0j . (C2) 



S'')(-^\J^^) 



Also, x^- (i)a;^ (t), which obtains from the interchange i <-> j, is the same. On the other hand, 

(1) (2) , (1) (2) 1/ (1) , . (l)w (2) . (2)n , , ,^o\ 

y> 'y)' + z\ 'z)' ^ -{y\ + iz\ >){y) > -iz)')+ h.c. (C3) 

To calculate this noise average, we make use of Eas. (|67l72|l . and 



{T)xf\T') = V^COS</),;e-5(--')0(^_^')^ 



0^{T)yf\T') = -2V^sin(/.,e-^("'-")0(T-T')cos2cc;(T-T'). (C4) 

where the stochastic noises Oi are defined as dWi = 6idt, and 9{t) is the step function. Performing the necessary time 
integrations and discarding 0{l/uj) terms, we get 

4 



(1) (2) , (1) (2) 'i . 2 J, • 2 

A y) + zl 'z) ' = -r]irij sm 4>i sm 
Finally, the average relative purity between the two single observer density matrices reads 



Vi Vj + zl 'z)j = -r]irij sin 4>i sin 0^. (C5) 



Otj = 2 + ViVj 



cos^ (pi cos^ (/)j H — sin^ (f>i sin^ (/)j 
9 



(C6) 
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